INCOMPRESSIBLE LIMIT FOR THE COMPRESSIBLE FLOW OF 

LIQUID CRYSTALS 



DEHUA WANG AND CHENG YU 

Abstract. The connection between the compressible flow of flquid crystals with low 
Mach number and the incompressible flow of liquid crystals is studied in a bounded 
domain. In particular, the convergence of weak solutions of the compressible flow of 
liquid crystals to the weak solutions of the incompressible flow of liquid crystals is proved 
when the Mach number approaches zero; that is, the incompressible limit is justified for 
weak solutions in a bounded domain. 



1. Introduction 

In this paper, we consider the incompressible limit of the following hydrodynamic system 
of partial differential equations for the three-dimensional compressible flow of nematic 
liquid crystals [21 [HI [28]: 

Pt + div(/)u) = 0, (1.1a) 
{pu)t + div(pu 5) + VP{p) = JlAu - Adiv (vd Vd - (^|Vdp + F(d))/3^ , (1.1b) 

dt + 5- Vd = ^(Ad-/(d)), (1.1c) 

where p > denotes the density, 5 E the velocity, d € the direction field for the 
averaged macroscopic molecular orientations, and P = op is the pressure with constants 
a > and 7 > 1. The positive constants Jj,,X,0 denote the viscosity, the competition 
between kinetic energy and potential energy, and the microscopic elastic relation time for 
the molecular orientation field, respectively. The symbol denotes the Kronecker tensor 
product, I3 is the 3x3 identity matrix, and Vd Vd denotes the 3x3 matrix whose 
ij-th entry is < dx^djOx^d >. Indeed, 

Vd Vd = (Vd)^Vd, 

where (Vd)"'' denotes the transpose of the 3x3 matrix Vd. The vector-valued smooth 
function /(d) denotes the penalty function and has the following form: 

/(d) = V^F{d), 

where the scalar function F{d) is the bulk part of the elastic energy. A typical example 
is to choose F(d) as the Ginzburg-Landau penalization thus yielding the penalty function 
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/(d) as: 



F(d) = -L(id|2 - 1)^ /(d) = (id|2 - i)d, 

40-0 ZCJq 



where ctq > is a constant. We refer the readers to [U EJ [l2l [T71 [261 EE] for more physical 
background and discussion of Uquid crystals and mathematical models. 
Recall that the Mach number for the compressible flow (II. ip is defined as: 



P'{p) 

Thus, letting M approach to zero, we hope that p, d keep a typical size 1, u of the order 
e, where e G (0, 1) is a small parameter. We scale p, u, and d in the following way: 

p = Pe{et,x), u = eue{et,x), d = de{et,x), 

and we take the viscosity coefficients as: 

Jl = efie, X = e'^Xs, = e6e, 

where the normalized coefficients /ig, A^, and 0^ satisfy 

fig ^ fx, — )■ A, 6e ^ 9 as e — 0"*", 

with ^, A and 9 positive constants. Under this scaling, system (jl.ip becomes 

^ + div(/5,u,) = 0, (1.2a) 
^%^+div(p,u,0u,)+vl/,J 

Ot E'^ 

= fle^U, - A,div (vd, Vd, - (^|Vd,|2 + F(d,))/3) , (1.2b) 

f)d 

^ + u,-Vd, = 0,(Ad,-/(d,)), (1.2c) 

where we take a = 1 because the exact value of a does not play a role in our paper. The 
existence of global weak solutions to (jl.2p in bounded domains was established in \i2\ I31j. 
By the initial energy bound (j2.10p below, we can assume that the initial datum is of 
the order 1 + 0(e), so it is reasonable to expect that, as e — )■ 0, — 1 and (|1.2ap yields 
the limit divu = 0, which is the incompressible condition of a fluid, and the first two terms 
in ()1.2bp become 

Ut + div(u u) = ut + (u • V)u. 
The corresponding incompressible equations of liquid crystals are: 

Ut + u • Vu + Vvr = ^Au - Adiv(Vd Vd), (1.3a) 
dt + u- Vd = e(Ad-/(d)), (1.3b) 
divu = 0. (1.3c) 

Thus, roughly speaking, it is also reasonable to expect from the mathematical point of 
view that the weak solutions to ()1.2p converge in suitable functional spaces to the weak 
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solutions of (jl.3p as e — )• 0, and the hydrostatic pressure vr in (jl.3ap is the "hmit" of 



in (|1.2b|) . This paper is devoted to the rigorous justification of the convergence of the 
above incompressible hmit (i.e., the low Mach number limit) for global weak solutions of 
the compressible equations of liquid crystals in bounded smooth domains. We remark that 
the existence of global weak solutions to the incompressible flow of liquid crystals ()1.3p 
was established in Lin-Liu 1291. 



When the direction field d does not appear, (II. 2p reduces to the compressible Navier- 
Stokes equations. Lions-Masmoudi |34] investigated the incompressible limits of the com- 
pressible isentropic Navier-Stokes equations in the whole space and periodic domains using 
the group method generated by the wave operator, a method introduced in earlier works 
[39] which requires certain smoothness of solutions. The study in bounded smooth 
domains with the no-slip boundary condition on the velocity is much harder than that in 
the whole space or periodic domains, because in bounded domains, there are extra dif- 
ficulties arising from the appearance of the boundary layers, and the subtle interactions 
between dissipative effects and wave propagation near the boundary, and hence requires 
a different approach. Desjardins-Grenier-Lions-Masmoudi in [6] relied on spectral analy- 
sis and Duhamel's principle to treat these difficulties to find the limit of global solutions 
in a bounded domain. These results have been extended by others; see for examples 
[21 Uni O |35l m] . We also remark that in Hoff [20] some convergence results were proved 
for well-prepared data as long as the solution of incompressible limit is suitably smooth. 
For the case of nonisentropic flows, see [13\ [T^ for some recent developments. Recently, 
Hu-Wang [22] studied the convergence of weak solutions of the compressible magnetohy- 
drodynamic equations to the weak solutions of the incompressible magnetohydrodynamic 
equations when the Mach number goes to zero in a periodic domain, the whole space, 
or a bounded domain; and Jiang-Ju-Li [23] studied the incompressible MHD limit in the 
inviscid case as long as the strong solution of the incompressible inviscid MHD exists with 
periodic boundary conditions. For other related studies on the incompressible limits of 
viscous and inviscid flows, see [B [1 [III [II [2ll [Ml IIS [271 [36l ISZl [Ml 112 and the references 
in |13j . Finally, we remark that the incompressible flow can also be derived from the 
vanishing Debye length type limit of a compressible flow with a Poisson damping; see for 
examples ^ |8j . 

Li this paper, we shall establish the incompressible limit of ()1.2p in a sufficiently smooth 
bounded domain C M^. As mentioned earlier, this limit problem in a bounded smooth 
domain has more difficulties and requires a different approach due to the appearance 
of the boundary layers and the subtle interactions between dissipative effects and wave 
propagation near the boundary. Comparing with those works on the compressible Navier- 
Stokes equations, we will encounter extra difficulties in studying the compressible liquid 
crystals. More precisely, besides the difficulties from compressible Navier-Stokes equations, 
the appearance of the direction field and the coupling effect between the hodrodynamic 
equations and the direction field should also be taken into account with new estimates. 
We will overcome all these difficulties by adapting the spectral analysis of the semigroup 
generated by the dissipative wave operator, Duhamel's principle, and the weak convergence 
method to establish the convergence of the global weak solutions of compressible flow of 
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liquid crystals (|1.2p to the weak solutions of the incompressible flow of liquid crystals (jl.3p 
as e goes to zero in a bounded domain. 

We organize the rest of the paper as follows. In Section 2, we will provide some prelim- 
inaries and state our main result. In Section 3, we will prove in four steps the convergence 
of the incompressible limit in a bounded domain. 



2. Preliminaries and Main Results 

We consider the incompressible limit in a smooth bounded domain C M'^. To state 
precisely our main result, we need to introduce a geometrical condition on n (cf. ^). Let 
us consider the following overdetermined problem: 

— Aip = \ip in fi, — — = on and w is constant on dVt. (2-1) 

ov 

A solution to ()2.ip is said to be trivial if A = and is a constant. We say that Vt satisfies 
the assumption (H) if all solutions of ()2.ip are trivial. In the two-dimensional case, it was 
proved that every bounded simply connected open set with Lipschitz boundary satisfies 
(H). We refer the readers to [6] for more information about assumption (H). 

Let us recall the definition of Leray's projectors: P onto the space of divergence-free 
vector fields and Q onto the space of gradients, defined by 

u = Pu + Qu, with div(Pu) = 0, curl(Qu) = 0, (2.2) 

for u € LP'. Indeed, in view of the results in jl5], we know that the operators P and Q are 
linear bounded operators in W^'^ for s > and 1 < p < oo in any bounded domain with 
smooth boundary. 

We consider a sequence of weak solutions {(pe, u^, d£)}e>o to (|1.2p in a smooth bounded 
domain C with the following boundary condition: 

Ue|9o = 0, d,|ao = d°, (2.3) 



and initial condition: 



satisfying 



Pe\t=0 = Pe, Pe'as\t=0 = m°, de\t=Q = d^, (2.4) 



>0, pUL-'m, (2.5) 



27 



G L~ (n), = if = 0, (2.6) 

pO|u°l'eLi(J]), d'^.GHHn), (2.7) 

converges weakly in to some uq as e — ?■ 0, (2-8) 

d^ converges weakly in to some dg as e — )• 0, (2-9) 



and 



/J^P°|u°P + ^A,|VdOp + A,F(dO; 



dx 



(2.10) 



+ ^2(7^1) /„ i(p'r - IP's + (7 - 1)) dx < c, 
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for some constant C > 0. We remark that (j2.10p implies, roughly speaking, that is of 
order 1 + 0(e) since 

{plV - IpI + (7 - 1) = iplV - 1 - 7 (P° - 1) , 
and p"^ is a convex function for 7 > 1. As proved in \42\ I31j . for any fixed e > 0, there 
exists a global weak solution (/Je,U£,d£) to the compressible flow of liquid crystals (jl.2p 
satisfying 

p,GL-([0,r];L^(l})), 
Vp7u, GL~([0,r];L2(J7)), 
u, GL2([0,r];Fi(f])), 
& L\%T]-H\^))rM^{[Q,T]-H\^)), 
for any given T > 0; and in addition, 

p,u, GC([0,r];L^(O)), 
p,GC([0,r];LP^(f])), 

if 1 < p < 7; as well as 

Ee{t)+ps[ [ \Vu,\^dxdt + Xeee [ [ {Ad^ - f {(1^)1^ dxdt < £,(0), (2.11) 

Jo in Jo Jn 

for t G [0, T] a.e., where 



and 



We now recall the existence result of global weak solutions to the incompressible flow 
of liquid crystals in [29] : 



Proposition 2.1. For uq G L'^{Q) and do G H'^{^}) with dolao G H^/'^{d9), system ((L3 
mi/i ^/le following initial and boundary conditions: 

u|i=o = uo(x) mi/i divuo = 0, d\t=o = do(x), 

and 

ubf7 = 0, d|af7 = dQ{x) 
has a global weak solution (u, d) such that 

u G L^{Q,T-H^{Q))r\L°°{^,T-L^{Q)), 

d G L2(0,r;i/2(f^))nL°°(0,r;i/i(S7)), 

and 

- /" (|u|2 + A|Vd|2 + 2AF(d))(i2;+ ! [ {plVul"^ + Xe\Ad - f{d)\'^) dxdt 
2 Jfi Jo Jn 



<\ I (|uo|2 + A|VdoP + 2AF(do))dx 
^ Jn 

for all T G (0,oo). 
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Remark 2.1. The global weak solutions obtained in Proposition l3.1l are the weak solutions 
in Leray's sense. The uniqueness of such solutions can be reached in two-dimensional 
spaces, see for example [29J. For more details on the existence and regularity of weak 
solutions to the incompressible flow of liquid crystals, we refer the readers to \29\ [30l 118] , 

Our main result reads as follows: 

Theorem 2.1. Assume that {(p^, Ug, d£)}£>o is a sequence of weak solutions to the com- 
pressible flow of liquid crystals (jl.2p in a smooth bounded domain Q CM.^ with the initial 
and boundary conditions (j2.3p - (|2.10p and 7 > |. Then, for any given T > 0, as e ^ 0, 
{(/>£, Ue, de)} converges to a weak solution (u, d) of the incompressible flow of liquid crys- 
tals (|1.3p with the initial data: u.\t=o = Puq, d\t=o = do and the boundary condition: 
u|ar2 = 0, d\QQ = do. More precisely, as e — 0, 

Pe converges to 1 in C{[0,T]; L'^ (n)); 

U£ converges to u weakly in L^((0,T) x ft) and strongly if VL satisfies condition (H); 

de converges to d strongly in L^{[0,T]; (ft)) and weakly in L'^{[0,T]; (ft)). 



3. Proof of Theorem 12.11 

In this section we prove Theorem 12.11 in four steps. 

3.1. A priori estimates and consequences. We first recall the spectral analysis of the 
semigroup [6j generated by the dissipative wave operator. Let {A^Qj^gN {^k,o > 0) be 
the nondecreasing sequence of eigenvalues and {$fc,o}fcGN in L'^{^) be the eigenvectors with 
zero mean value of the Laplace operator satisfying the homogeneous Neumann boundary 
condition: 

-A^k,o = A|o^fc,o in ^, = on dQ, 

where v is the unit outer normal of $7. By the Gram-Schmidt orthogonalization method, 
it is possible to assume that {^fc.olfceN is an orthonormal basis of and that up to a 

slight modification, if Afc,o = M,o-, and k I, then 



V$fc,o • V$;,orfs = 0. 

I an 

From (|2.1ip and the conservation of mass, we have for almost all t > 0, 

^ Qft|u,|2 + ^2^^^_i) (p2 - 7Pe + 7 - 1) + Y|V,d,|2 ^ 
+ He [ [ \VMe\^dxdt + \eee I [ \Ade - fide)\'^dxdt 

Jo Jn Jo Jn 

I (^P°|u°P + i2(:^((/'°)" - 7/^° + 7 - 1) + ||V,d,T + AF(d°)) dx 



(3.1) 



< 

< c. 

By (j3.ip . we have the following properties: 

is bounded in L'^{[0,T]; L'^{n)), (3.2) 
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d, is bounded in L^{[0,T]; H\Q)), (3.3) 
^ -(pj -7P£ + 7-1) is bounded in L°°([0,r];Li(O)); (3.4) 



£2(^-1) 

and 

Vu. is bounded in L^{[0,T]; L^{n)), (3.5) 
Ade - f{de) is bounded in L^{[0,T]; L^{n)) (3.6) 

for all T > 0. 

Let us recall the following basic fact: for some cq > and for all x > 0, 
x"^ - 1 - 7(x - 1) > co\x - 1|^ if 7 > 2, 
x*^ — 1 — 7(x — l)>co|x — l|^if7<2 and x < R, 
x"^ — 1 — 7(x — 1) > co|x — ip if 7 < 2 and x > R, 
where R G (0, cxd). Thus from (j3.4p . we have 



(y^\P^ - lPx[p,-l|<l/2 + ^\pe - irXlp,-l|>l/2) dx < C, 



(3.8) 



where x is the characteristics function and C denotes a generic positive constant hereafter. 
By (j3.8p . one has 

sup \\pe - l||L7(m < Ce""/^, and sup \\pe - l||L«(m < C'e, 

t>0 t>0 

where k = min{2,7}, which implies that 

Pe^lmC{[0,T];L^{Q)) as e ^ 0. (3.9) 
Now we split the velocity 

= + with = U£X|p,_i|< 1 , = u^X|p,_i|> 1 , 
which means that 

sup / |u]|2dx<2sup / pe\ue\'^dx <C (3.10) 

and 

\\^e\\h{n) ^ jjP^ ~ lll"^!^'^^ ^ (^e||Ue||^2./(._i)(f^) < Ce\\VUe\\l2^^), (3.11) 

where we used some embedding inequality. By (j3.10p and (j3.1ip . it is easy to see that 

is bounded in L^{[0,T]; L^{n)), 

and 

u^e-i/^ is bounded in L'^{[0,T]; L^{n)). 

This implies that is bounded in L'^({0,T) x $7) for all T > 0. 
By ()1.2ap and ()3.9p . one deduces that, 

divu^ weakly in L'^{[0,T];L'^{n)) 

for all T > 0, where we used a fact 



aldivuepdxdt <C [ I \Vus\'^dxdt < C. 
_ _! Jo Jn 
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By (j3.6p . smoothness of /, and the standard ehiptic theory, we have 

V'd,€L\[0,T];L\n)). 

On the other hand, multiplying on the both sides of (jl.2cp . then applying maximal 
principle, we have G L°°{[0,T] x Q). Thus, 

d,GL\[0,T]-H\n)). 

Similarly to [42j, using the Gagliardo-Nirenberg inequality, 

1 i_ 

l|Vds||i4(n) < c||Ad^|||2(f^)||d^|||^(^) +c||d^||ioo(f^), 

one has 

Vde G L^((0,r) X 0). (3.12) 
Summing up the above estimates, we can assume that, up to a subsequence if necessary, 

Us-^u weakly in L'^{[0,T];H^{n), 

divue weakly in L'^{[0,T]; L'^{n)), 
in C([0,r];L^(J7)), 
Ad, - /(d,) ^ Ad - /(d) weakly in L\[0,T]; L\n)), 

and 

de^d weakly in L°°{0,T; H\n)) D L^{[0,T]; {fl)). 

To show the strong convergence of d,, we rely on the following Aubin-Lions compactness 
lemma (see [53]): 

Lemma 3.1. Let Xq,X and Xi be three Banach spaces with Xq Q X Q Xi and Xq and 
Xi be reflexive. Suppose that Xq is compactly embedded in X and that X is continuously 
embedded in Xi. For I < p,q < oo, let 

W=l^ue LP{[0,T];Xo) : ^ € L^([0, T]; Xi)| . 

Then the embedding ofW into ^^([0, r];X) is also compact. 
From ()1.2cp . it is easy to see that 

\\dtds\\L2(n) < WOsi^ds - f{ds))\\L2{n) + \\ue ■ VdsWi^n) 

< C||ue||i4(f,) + C\\Vd,\\l,^^^ + C7||Ad, - /(d,)|U2(f,), (3.13) 

< C||Vu,||i.(^) + C||Vd,||i4(^) + C||Ad, - /(d,)||i.(f,). 
This, combined with (|3.5p . ()3.6p . (j3.12p and embedding theorem, implies that 

\\'^tde\\L2{[0,T];L'^(n)) < C. 

Since C C and the injection M> is compact, we apply Lemma 13.11 to 
deduce that the sequence {d,} is precompact in L^(0, T; H^{Q)). By taking a subsequence 
if necessary, we can assume that, 

d, ^ d weakly in L'^{0,T; H'^{n)), 
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aoad 

^ d strongly in L'^{0,T- H\Q)). 

Therefore, by a standard argument, we deduce that the Hmit u and d satisfy equation 
()1.3cp in the sense of distributions. By strong convergence of d^ and smoothness of F, we 
deduce the convergence of the nonhnear term 

Vd,0Vde-(^|Vd,|2 + F(d,))I 

^ Vd0 Vd- (i|Vd|2 + F(d))/ ase^O 
in the sense of distributions. Denoting 

vr. = ^pJ-y|Vd,|2-A,F(d,), 
then we rewrite equation ()1.2bp as 

^^g"^^ + div(p,u^ Ue) + VvTe = /i^Au^ - Aediv(Vde Vd^). (3.14) 
We project equation (j3.14p onto divergence-free vector fields: 

dtP{peU,) + P{dw{peU, U,)) - fleAPUe = -P(Aediv(Vd, Vd,)), (3.15) 

where P is defined by (12. 2p . Then (I3.15P yields a bound on dtP{pe'yie) in 

L'{[o,T]-H-\n)) + L\[o,Ty,w~''Hn)) + L\[o,T]-H-'{n)), 

and hence in {[0,T];W^^'^ (Q,)). In addition, P{psU^) is bounded in 

L'^i[0,T];L^in))nL\[0,T];Vin)) with - = - + ^. 

r 7 D 

To continue our proof, we need the following lemma (cf. Lemma 5.1 in [32]). 

Lemma 3.2. Let the functions gn, hn converge weakly to the functions g, h, respectively, 
in LP^{0,T;LP^{n)), L^i (0, T; ^^2(17)), where 1 < pi,p2 < oo and 

1 1 _ 1 1 _ ^ 
Pi qi P2 q2 



Assume, in addition, that 

is hounded in L^[Q,T;W~'^'^{^)) for some m>0 independent of n, 

and 

\\hn - hn{t, ■ + C)llL9i(o,r;L92) — > fls |C| — > Uniformly in n. 
Then gnhn converges to gh in the sense of distributions in (0, T) x Q. 

Applying Lemma [3T2l we deduce that P(p£Ue) • i-'(ue) converges to |up in the sense of 
distributions. It is easy to see that Pu^ — )• u = Pu in the sense of distributions because 
the weak convergence of P\i^ to u = Pu in L^((0, T); L^(0)) and 

/ / i\P^e\'^ - Pipers) ■Pu,)dxdt < C\\pe - 1 lb([0,T];LT(n)) ||Ue ||i2([o,T];L»(n)) ' 
"/ 

with s = < 6 since 7 > f • 
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3.2. The convergence of Qu^. To prove our main result, it remains to show the conver- 
gence of the gradient part of the velocity Qu^. The argument for the convergence of Qu^ 
in our paper follows the same line in |6] and |22j . except the argument for the direction 
field dg. For the convenience of readers and the completeness of argument, we provide the 
details here. 

First, we introduce the spectral problem associated with the viscous wave operator 
in terms of eigenvalues and eigenvectors of the invisicid wave operator L. In the sequel, 
we write the density fluctuation as 



1 



e 



and 



and 



($,m) . We define the wave operators L and in x Viyt) as follows: 

m 



L 



divm 
V<I> 



m 



L 



m 



+ e 







The eigenvalues and eigenvectors of L read as follows: 



fc,0 



m 



fc,0 



. V'^'fc.o 



(3.16) 
(3.17) 

(3.18) 



and 



Lcfy^Q = ±iXkfl(j)f Q in 17, m^Q • i/ = on dQ. 
And we need the following lemma [6]: 

Lemma 3.3. Let Q, he a hounded domain ofW^ and let k > 1, N > 0. Then, there 
exist approximate eigenvalues i^^^j^ and eigenvectors (jy^^j^ = e ad "ifc ^ at)^ o/ 
such that 



''''^k,e,N'^k,e,N + ^k,e,N^ 



with 



iX 



k,e,N 



±iXk,o + i>^ki^ + 0{e), where Re{iX^^) < 0, 



and for all p G [l,oo], we have, 

< Cp{^ef+^'P and cPk,e,N - 



^k,e,N 



LP(n) 



LP{n) 



i/p 



Remark 3.1. The key idea is to construct an approximation scheme of in terms of 
(fy^Q- We refer the readers to [6] for more details and the proof. From Lemma 13.31 and its 
proof, we have, for any integers k, 



which satisfies 



Re{iX^-^) < 0. (3.19) 
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We observe that the first order term iX^i clearly yields an instantaneous damping of 
acoustic waves, as soon as Re{iX^^) < 0. Thus, we let 

/ C N 

to be a collection of the all eigenvectors ^k,o of the Laplace operator such that 

Re{iXf^) < 0. 



Denote 

that is to say, when k £ J, we have 



J = n\i, 

Re{iX^^) = 



due to (|3.19|) . This implies A^^o = 0- Iii the case that A^^o = 0; "^fe o ™^st vanish on dQ 
and therefore not only m^Q-z/ = but also = on d^l. Thus, no significant boundary 
layer is created, and there is no enhanced dissipation of energy in these layers. 

Remark that { jfcgN is an orthonormal basis of L^{Q) functions with zero mean 
value on i}. We write 



fcGN 

where 



Afc,0 



if,g) = / f{x)g{x)dx. 
Jn 

We split Qug into two parts Qiu^ and Q2Ue, defined by 

V At. n J Xl 

and 



A 



k,0 



keJ 

which, respectively, correspond to damped terms and nondamped terms. What remains 
is to show that 

QiUs^O in l2([0,T] X 

and 

curldiv(Q2in£ ^ Q2Ue) —?- 

in the sense of distributions if J 7^ 0, which is equivalent to saying that div{Q2m^ Q2^e) 
converges to a gradient in the sense of distributions. 

Let us observe that in view of the bound on Ug in L^(0,T; Hq{Q,)), the problem reduces 
to a finite number of modes. Remark that the eigenvalues {A|o}fc>o is a nondecreasing 
sequence, we have 
rT 

E 



k>N 



C 

dt < -2 |Vu^|i2((o,T)xn) ' i = 1 or 2. 

^Af+l 
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Letting — )• oo, then — )• oo, which imphes that {Qiu^,m^Q) — )• in L^(0, T) for any 
k> N. So we need to show that (QiU£,m^Q) converges to strongly in L^(0,T) for any 
fixed k and study the interaction of a finite number of terms in dw{Q2U^ ® Q2^e)- 
RecaUing that ipe = , we have 

and 



e\{Q{ipeUe),V^ 



k,0) 



(3.20) 



which tends to zero in L^(0, T) due to 7 > |. So we move to study {Qm^, iii^q). We write 



with 



:{t) = ( 1^ ) . (3.21) 



It is easy to see 

2(Qm„mJo)=/3±,-/3,^,. 
Applying Lemma 13.31 with N = 2 and the Holder inequality, one obtains 



(3.22) 

<Ce2 ( ||(y9e||Loo(ro,Tl;L'=(C)) + llniell 

where 

11 1 , 

a = min] 1 , 1. 

^ k'2 27-^ 

It remains to show that bf^{t) = {4>e{t)i'Pk eo) converges to zero strongly in L^([0, T]) 
when k £ I and check the oscillations when k G J. 

Using L* to denote the adjoint operator of with respect to (•,•), we have 

where 

ge = -div(me (g) rrie) - Vvr^ - A£div(Vde Vd^). 
Taking the scalar product of (j3.23p with (p^^ 2) one obtains 

idt<t>e, c/>t,2) - '^'^,2) = <eit) (3-24) 

where 4^^{t) = (5^, m^^^^g) + ^^,2)- 
Letting b^^{t) = {(ps, 4>^^ 2) observing that 



INCOMPRESSIBLE LIMIT OF LIQUID CRYSTALS 



13 



we rewrite (j3.24p as 



1- 



-iXt..bt. = ct^. (3.25) 

We can view p.25p as an ordinary differential equation, which satisfies the initial value: 

b^Jt=o = 6j,(0). (3.26) 

3.3. The case k & I. By the basic theory of ordinary differential equation, the solution 
to (|3:25l) - (l3:26|) is given by 

ft 



iX 



iX 



(3.27) 







Next, we need to estimate it). Note that 



which is helpful in estimating the first part of (I3.27P : 



Wbl.m'^'-'-'nL^io^T) < C\\bl^me'^''^''"^\\LHo,T) < Cet, 

where we used Lemma 13.31 

To estimate the second term in (j3.27p . observe that for any a £ L'?(0, T) and I < p,q < 



± N t 



(3.28) 



oo such that - + - = 1, we have 



p q 



f 

Jo 



JX 



k,£,2 - 



'a{s)ds 







Jo 





\a{s)\ds < C\a\Li{o,T)V£'' ■ 



(3.29) 



To prove 



6±^(t) ^ strongly in L\[0,T]), 



it remains to show that c^^ is bounded in L'?(0,T) for some q > 1. Remark that 



where 



ci 



< Cl + C2 + C3 + C4, 



/ (me (g) ms){t) ■ Vm^_, ^dx 



and C2,C3,C4 will be defined as follows. Using 



and integration by parts, we have 

ci{t) < e l^PeUe Us ■ Vm^^ ^\dx + C ju^ • Vu^ • m^^ ^\dx 
Jn ' ' Jn 

I 111 2ii 

- ll™fc,e,2lli°°(n)l|Ue + UellL2{n)l|VUe||L2(C) 

+ £||¥'e||L°°([0,T];L-(f7))||(Ue)^||i«./«.-i(f^)||||Vm^^^^2lli°°(n) 

< C||Ue||L^([o,T];L2(n)) • ||Vue||i2(f^) + Ce^^^ || Vu^ || ^2(q) + Ce^/^||Vu£||i2(f^). 
For the second term 

C2 = / (TTedivm^^ 2)dx, 
Jn 



(3.30) 
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recalling that 
we have 

C2(i) < ||7I"e||L°°([0,T];Ll(n))(ll^fc,e,2llL°°{n) + II -^^,£,2 II (f^) ) ^ ^- (3.31) 

The third term can be estimated as: 

^ £~"^ll-^fc,£,2llL-/(— i)(n)ll<^£llL°°([0,T];L«(C)) (3.32) 
< (^gl/2-l/2K_ 



And the last term is 



C4(t) 



/ div(Vd£ Vd£)m^^ 2*^^ 



<C JjVd,\'\Vui^^Jdx (3.33) 

< C||VmJ^^2lU°°{n) II Vd^ 11^2(0) 
<C7||Vd,|U2(o), 

where we used Vde G ^^(J^). From (I330]l - (l3:33]l . it is easy to see that the term c. is 
bounded, and consequently, 

^k,e ~^ ^ Strongly in L^(0,T). 
To complete our proof, we need to consider further the case k G J. 

3.4. The case k € J. As in Remark 13. H when k £ J, we have A^^ = 0. This implies 
that, together with (|3.'27p . 

g±iAfc,oi/£5±^ is bounded in L'^{0,T), 

and 

Q^(g±iAfe,oV£5±j is bounded in ^/lL^{0,T) + 1^(0, T) for some p>l. 

It follows that, up to a subsequence if necessary, b^^ converges strongly in L^(0, T) to 
some element bf^^^. Since /Og — >• 1 in C{[0,T]; L"'{i})) and b^^ is uniformly bounded in 
L^(0, T), one obtains that 



in the sense of distributions. Thus, we only need to consider the terms 

btsKs^^^ for all A:,/ G J. 
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On the other hand, 



Afc,0 A; 







in the sense of distributions, where we used the strong convergence of 

g±Afc,ot/£5±^ in L2(0, T) when A; G J. 

If ^kfi = Ai^O) we can write 

div(V^>fc,o ® V$/,o + V«>,,o ® V«>fc,o) = -A|oV($fc,o^z,o) + V(V$fc,o • V«>,,o), 
which is a gradient, and thus disappears in the pressure term. It remains to consider the 
case Aa:,o / A;^o- In this case, we have a fact that 

^tosMbtscit) ^ L\[0,T]) 

for all k,l £ J due to 

btscit) ^ L\[0,T]) for ah A: G J. 
This implies, together with the Riemann-Lebesgue Lemma, 

r g^(A.,o-A,o)Veb±^^^(i)5± ^^(i)rfi ^0 as e ^ 0. 
Jn 

Thus, we conclude that 

in the sense of distributions. Hence, the finite sum, as k,l < N, 



converges to a gradient in the sense of distributions, which means that 

div {psQ2Ue (gi Q2Ue) 

converges to a gradient in the sense of distributions. 
The proof of Theorem 12.11 is now complete. 
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